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Abstract 



We show that the sensitivity of wave speckle patterns in disordered non- 
linear media to changes of scattering potential increases with sample size. 
For large enough sample size this quantity diverges, which implies that at 
given coherent wave incident on a sample there are multiple solutions for the 
spatial distribution of the wave's density. The number of solutions increases 

exponentially with the sample size. 
Suggested PACS index category: 05.20-y, 82.20-w 
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If a coherent wave described by a field <p(r, e) propagates in an elastically scattering 
medium, the spatial dependence of its "density" n(r,e) = \(p(r,e)\ 2 exhibits speckle: n(r,e) 
is a random, sample specific function of coordinate r. Here e is the wave's energy. In the 
cases of noninteracting electrons and electromagnetic waves propagating in linear media the 
theory of sensitivity of speckle patterns to a change in scattering potential was developed 
long ago ' 1_5 l It was shown that the sensitivity is very large, but finite. 

In this article we consider the same question in the case where a wave propagates in non- 
linear media. For the sake of concreteness we consider the situation where the propagation 
of the wave is described by a nonlinear Schrodinger equation 



Here m is the wave mass, u(r) = (3n{r) is the effective nonlinear potential and u(r) is 
a scattering potential which is a random function of the coordinates. Similar equations 



of hydrodynamic turbulence ^ , and the theory of turbulent plasma t 8 ' . We will assume white 
noise statistics in u(r): (u(r)) = , (u(r)«(ri)) = -rh^r — r-i). Here brackets () correspond 



Let us consider the case where a coherent wave <f>o(r) = v /noexp(zfcr) with momentum k 
is incident on a disordered sample of the thickness L 3> I (See the insert in Fig.l). We 
will show that the sensitivity of the nonlinear speckle pattern n{r) to a small change in 
u(r) increases with sample size L. At arbitrarily small n and for arbitrary sign of /? the 
sensitivity become infinite provided L is large enough. This implies that at given coherent 
wave incident on a sample Eq.l has many solutions. This is very different from the case 
of uniform nonlinear media, where types of instabilities depend on the sign of (3. (See, for 
example, ^.) 

The r dependence of the average density (n(r,e)) can be described by the diffusion 
equation, which is equivalent to calculation of the diagrams shown in Fig. 2. a. We use the 
usual diagram technique for averaging over realizations of random potential ' 9 l If \(3n \ <C 




(1) 



appear in the theory of electromagnetic waves propagating in nonlinear media ' 6 1, the theory 



to averaging over realizations of u(r) and I is the elastic mean free path (I >> k 



(2em) - 5). 
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'f^ one can neglect the nonlinear corrections to the diffusion coefficient D — j^. We 
obtained this criterion by calculating the transport scattering cross section on the effective 
potential (3n{r). To do so we used calculated in I 11,4 ! spatial correlation functions of the 
density l 10 l In the case of the sample geometry shown in the insert of Fig.l, we have 
(n(r)) = n . 

We can characterize the speckle pattern n(r) and its sensitivity to a small change of scat- 
tering potential A-u(r) = u'(r)—u(r) by correlation functions {Sn(r)Sn(ri)) and K(r, ri) = 
((An(r)An(ri))). Here 5n(r) = n(r) — (n(r)); An(r) = n'(r,{u'(r)}) — n(r,{u(r)}), 
n(e, r, {-u(r)}) and n'(e, r, {li'(r)}) are solutions of Eq.l with scattering potentials -u(r) and 
u'(r) respectively, and the brackets (()) correspond to averaging over both realizations of 
u{r) and realizations of Au{r). We will assume that ((u(r)u' (ri))) = £/ 2 exp(^p^). To 
get the value of K(r, j*i) at \r — r±\ ^> I and in the first in A«(r) approximation one can 
generalize the Langevin approach for calculations of mesoscopic fluctuations t 4 ' 11 ' to include 
nonlinear effects 

^-5J(r)=0; 5J(r) = -D^-dn(r) + J ext (r); (2) 

J ex t(r, Mr)}) = J ext (r, {«(r)» + / dr' 5 -j^- W) + Au(r')) (3) 

2irl 

(JUrUUri)) = ^(n(r)) 2 5(r - Pl )* y (4) 

<ww > = ^ (r - ri){G( ^ yi)Wr)>x 

((n(r;))G(r',r) + (n(r')>G(rl, r)) - (n(r'))(n(r' 1 ))G(r',r))G(r' 1 ,r)} (5) 



Au(r) = u'(r) - u(r) = (3An(r) (6) 
where G(r, ri) is the Green function of the equation 

-^G(r,r 1 ) = 5(r-r l ); (7) 



J ( r ) = 2k /m ^*( r )^^( r ) is the current density, 5J{r) = J(r)-(J(r)), J ext (r , {u(r)}) 
is a random external current source, (J ex t(r)) = ( ^ff ) = (JLti 7 ") 5 '^^ ) = °> and 
«, j are the coordinate's indices. Eqs.2,3,7 require the usual diffusion boundary conditions: 
G(r, r') = 0, < n(r) >= n at x = and n • ^G(r, r') = n • £{n(r)) = n ■ J ext (r) = at 
the closed sample's boundaries. Here n is a the unit vector normal to the boundary. 

Eqs.2-7 are a closed system which differ from ' 4,11 1 by the term in Eq.3 proportional to f3. 
They are equivalent to the summation of diagrams shown in Fig.2b-g. Diagrams, shown in 
Fig.2h, are responcible for the small nongaussian part of the distribution function of SJ & ^^ ■ 
They are proportional to a small parameter p^- <C 1 in the three dimentional case (d=3) 
and can be neglected. All diagrams responcible for localisation effects can be neglected as 
well. 

Let us first consider the linear case f3 — 0, u{r) = 0. Index (0) will indicate quantities 
calculated at (3 — 0. Solving Eqs.2-5,7 at \r — t*i| > I in d = 3 case we get M 

K^(r, ri ) = ((An^(r)An^( ri ))) ~ C~) 2 < Sn(r)Sn( ri ) >~ — ( 8 ) 

r/ /c z t|r — ri\ Tf 

where td = jj, and ^ = ^ characterizes the change in scattering potential. This can 
also be obtained by calculating the diagrams shown in Fig. 2b, c. The characteristic time 
~ "7J corres P on ds to a complete change in the speckle pattern due to the change of the 
scattering potential Att(r). One can get the same estimate from the requirement that an 



additional phase ~ , / L t m , which the traveling wave aquires due to the change in the 

V Dr f 

potential Au(r), is of order ir. If impurities have a cross-section of order p and they are 
shifted from their initial position by distances of order | then in the d = 3 case one has to 
change the positions of the = Llk 2 impurities to in order change the speckle pattern 
significantly ^ . Characteristic changes of energy Ae^* = and of the angle of incidence 
(see the insert in Fig.l), which change the speckle pattern significantly, can be 
obtained in a similar way l 1 ' 4 '. 

Let us now turn to the case (3 ^ 0. Expanding Eqs.2-6 up to second order in j3, and 
performing the average over realizations of u(r) and Au(r) in the d=3 case we get the 



correction to 7f(°)(r, ri). 

K^ir,^) ~ 1 K ( -°\r,r 1 ) (9) 

where 

7 = ( 3n^ )2( L )3 (1Q) 

The index (1) indicates quantities proportional to (3 2 . It can also be obtained by calculating 
the diagrams shown in Fig.2d-g or by estimating the additional phase which the wave trav- 
eling along a typical diffusion path will pick up due to the change in the effective potential 
A#(r) =pAn<®(r) 

«(Ax {1) ) 2 » = (^) 2 « / dsds 1 An^(r(s))An^(r(s 1 )))) ~ 7 «(x (0) ) 2 » (H) 

T 2 

Here integration is taken along typical diffusion paths of length j-. 

Eqs.9,11 imply that (((x {1) ) 2 >> > (((x {0) ) 2 )) and that Eq.l has many solutions. Let us 
estimate the number of the solutions in the D=3 case. It is convenient to expand 

u(r) = ■^='22m*u m n m (r) (12) 
yj L m 

over a complete set of eigenstates n m (r) of diffusion equation Eq.7 

d 2 

~ D ~^. nm ^ = E m n m{r) (13) 

where E m ~ r^ms are eigenvalues of Eq.13 and m = 1,2... labels the eigenmodes. Let us 
first regard u m as independent variables. The solution of Eq.l can be written as n{r) = 
n(r,Ui, ..Uk-.-)- Than using the selfconsistency equation u = (3n we get 
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7~2 m 3^ m = F m (u u ..u k ...) (14) 

where F m (ui, ...) = kL~ 1 riQ 1 m^l^ J drn(r,u\...)n m {r) are random sample specific functions. 

The problem of the investigation of properties of F m (ui, ...) as a function of u n is equiv- 
alent to the linear problem considered in I 1-5 !. To characterize the dimensionless functions 



F rn we calculate the following correlation functions with the help of Eqs.2-5,7: (a)mesoscopic 
fluctuations of modes with m ^ n are uncorrelated (SF m SF n ) = 0, where 5F m = F m — (F m ); 
(b) {{5F m f) ~ 1; and (c) 



([Frnjui, .-U n + Au n , ..) - F m (ui..U n , ..' 
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«fF„)»> - < As »> ' < 15 > 

Eq.15 means that the characteristic period of random oscillations of F m as a function of u n 
is of order unity, Au n ~ 1. In anticipation of these results we have introduced the factor 
mJ in Eq.12. 
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It is important that Eqs.14 with large enough m > M = 74 have unique solutions 

2 

u m ~ m~3 <g; 1, provided the values u m< u are given. We get the estimate for M from the 
requirement 7"2m3 <C 1. Therefore at 7 ~ 1 modes with m ~ 1 are the most important for 
determination of the number of solutions of Eq.l. This also follows, for example, from the 
long range r—r 1 dependence of K°(r, r±) and from the the fact that the main contribution to 
Eq.9 and the diagrams shown in Fig. 2 d-g is from integration over intermediate coordinates 
with \r — r'\ ~ L. Therefore we would like to introduce a model which captures the main 
features of the problem at 7 ~ 1: we put u m> \ = in the first of Eqs.14 for m — 1 and get 
the equation 

7 -ki = Fi(iZi,0,0,..) (16) 

It is equivalent to substitution in Eq.l u(r) — > ^■Uin 1 (r). (Then expanding Eq.l with 
respect to powers of (3 one can reproduce the values of the diagrams Fig.2b-g with the 
precision of the factor of order of unity.) In Fig.l we show a qualitative "graphical" solution 
of Eq.16 which corresponds to intersection of two functions: Fi 0, 0..) and 7~^Wi. It 
follows from Fig.l that at 7 > 1 both Eq.16 with a typical realization of the potential 
u(r), and consequently Eq.l, have many solutions. In this case the sensitivity defined as 
(^) 2 K(r,r') diverges. The main contribution to this divergency comes from realizations 
of u(r), when Fx(ux, 0, 0..) and are tangent to each other. The criterion 7 > 1 is 

equivalent to the inequality ((^sjj / n(r)dr) 2 ) > 1. In such a form this criterion is similar 
to the criterion of Stoner ferromagnetic instability in metals ' 13 '. 



We would like to mention that even at 7 < 1 there are rear realizations of u(r) which 
correspond to several solutions of Eq.l. Therefore, formally speaking, the sensitivity diverges 
at any 7. Obviously the conventional diagram technique is unable to describe the existence 
of many solutions of Eq.l. 

At 7 3> 1 the number of solutions of Eq.16 shown in Fig.l is of order 75. However, if 

7 ^> 1 not only ui, but also higher modes with 1 < m < M, are relevant. In this case 
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Eqs.14 have multiple solutions in the intervals \u m \ < ^m~3. Since both the amplitude of 
fluctuations and the periods in m — th direction of randomly rippled hypersurfaces F m (u\...) 
are of order unity, the number of solutions N of Eqs.14,1 is proportional to the volume of 

1 2 

the manifold \u m \ < 72 m^, m < M. As a result we have 

m M 2 3 

N ~ 1~ II m 3 = exp(a74) (17) 
1 

where a ~ 1. 

Similar phenomenon may occur in disordered metals with interacting electrons. The 
system can be unstable with respect to creation of random magnetic moments. In this case 
n(r) would play the role of magnetization density. This would correspond to Finkelshtain's 
scenario ' 12 l However, in this case to get a self consistency equation for n(r) we have to 
integrate over electron energies up to the Fermi energy, which decreases the amplitude of 
mesoscopic fluctuations of n(r). As a result, at small electron-electron interaction constant 
the situation with many solutions occurs only in the D=2 case and the characteristic spatial 
scale of integration over r will be of the order of the localization length in the linear problem. 
Thus the problem of interacting electrons in disordered metals remains unsolved. 

Above we considered the case when <j)(r,t) = <f>(r, e) exp(iet) is a complex quantity 
and n(r) = \<f)(r,t)\ 2 is time independent. Therefore the third harmonic, proportional to 
exp(3iet) is not generated. In the case of propagation of electromagnetic waves in nonlinear 
media 0(r) should be considered as a real quantity, which leads to generation of third 
harmonics. In this case the presented above consideration is valid only as long as the 
amplitude of the third harmonic is smaller than the amplitude of the first harmonic. It is 
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the case provided that <C 1. The letter criterion does not contradict to the requirement 
7 > 1. 

Finally, we would like to mention that the problem considered above is similar to the 
problem of classical chaos, where the the sensitivity of trajectories of motion to changes in 
boundary conditions exponentially increases with the sample size (See for example ' 14 ^). 
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FIG. 2. Solid lines correspond to Green functions of Eq.l with (3 = 0, dashed lines correspond 
to -r-jS(r — ri), the four solid lines vertices correspond to the factor /3, thick wavy line correspond 
to Au(r) and thin wavy lines correspond to density vertexes. 
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